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ABSTRACT 

I show here that the bulk of the dust mass in a galaxy can be equivalently estimated from: i) the full spectral energy distribution of dust 
emission, using the approach of Draine & Li (2007) that includes a distribution of dust grains and a range of interstellar radiation field 
intensities; ii) the emission in the wavelength range 100yun? < A < 500/im (covered by the Herschel Space Observatory), by fitting to 
the data a simpler single temperature modified blackbody. Recent claims on the contrary (Dale et al. 2012) should be interpreted as a 
caveat to use in the simpler fits an absorption cross section which is consistent both in the normalization and in the spectral index /3 
with that of the full dust model. I also show that the dust mass does not depend significantly on the choice of /?, if both the dust mass 
and the absorption cross section are derived with the same assumption on /?. 
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1. Introduction 



Oh. 



A full coverage of the Spectral Energy Distribution (SED) of 



"dust emission in galaxies has recently become available, thanks 
O mainly to the Spitzer Space Telescope (Werner et al. 2004) and 
£j to the Herschel Space Observatory (Pilbratt et al. 2010). Dust 
CZ) emission models can now extract from the infrared observations 
, ^ , a wealth of information on the dust composition, grain sizes and 
intensity of the dust heating sources. One such model is that 
i- — I of Draine & Li (2007, hereafter, DL07), which has been suc- 
J> cessful in reproducing both the global (Draine et al. 2007; Dale 
On et al. 2012, hereafter, D12) and the resolved SEDs (Aniano et al. 
ON 2012). 

A general result of the modelling is that emission for A > 
I lOOyum predominantly comes from dust heated at thermal equi- 
CN librium by a mean interstellar radiation field (ISRF); and that 
<0> this dust component constitutes the bulk of the dust mass in a 
CO galaxy (*98-99%; Draine et al. 2007, D12). If all dust grains 
\ I share the same size and composition, this emission is equiva- 
J> lent to that of a single temperature modified blackbody (MBB), 
? i.e. a blackbody multiplied by the dust absorption cross section. 
For a dust model including grains of different sizes and com- 
positions, and thus different absorption cross sections, the SED 
could be broader than that of a MBB, because different grains at- 
tain a range of thermal equilibrium temperatures. Thus, the mass 
obtained by fitting to the observed SED a MBB with an average 
absorption cross section could in principle be biased with respect 
to that derived with the DL07 approach, using the full dust grain 
model. Nevertheless Magrini et al. (2011), fitting Herschel data 
at A > 100/im for a sample of Virgo Cluster galaxies from the 
HeViCS programme (Davies et al. 2012), found that MBB dust 
masses are within »10% of DL07 masses, proving that the SED 
broadening due to the individual grain temperatures is minimal. 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Using the same spectral range for the SEDs of galaxies in 
the Herschel KINGFISH sample (Kennicutt et al. 201 1), D12 in- 
stead claimed that the MBB approach can underestimate the dust 
mass by up to a factor two, because "single blackbody curves do 
not capture the full range of dust temperatures inherent to any 
galaxy". 

Starting from the D12 dataset, in this note I derive indepen- 
dently the MBB dust masses. I confirm the results of Magrini 
et al. (2011) (Sect. 2). I also show that a meaningful deriva- 
tion of the dust mass can be obtained only if the same spectral 
index is used both in the MBB fitting and in the derivation of 
the absorption cross section (Sect. 3). Conclusions are drawn in 
Sect. 4. 



2. Dust masses from modified blackbody fits 

D12 measured the flux densities of the KINGFISH galaxies ob- 
served by Herschel with the PACS instrument at 70, 100 and 
160jt/m (Poglitsch et al. 2010), and the SPIRE instrument at 
250, 350 and 500/rni (Griffin et al. 2010). After complementing 
the Herschel observations with available infrared observations at 
shorter wavelengths, they derived various parameters character- 
izing the dust emission, including the dust mass, by using the 
DL07 approach. Here I use their photometry for the five Her- 
schel bands with A > 100/mi, and derive the MBB dust mass for 
the objects detected in all five bands (56 objects out of 61). The 
distances of the galaxies is taken from Kennicutt et al. (201 1). 

Under the assumption that all dust grains share a single tem- 
perature Td, and that the dust distribution is optically thin, the 
mass of dust Md can be estimated by fitting the observed flux 
densities f v to a MBB, 



fy 



Md 
D 2 



(1) 



where D is the distance of the object, B V (T&) the Planck func- 
tion, and K a bs the grain absorption cross section per unit mass (a 
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Fig. 1. Ratio between dust masses derived from MBB fits to the 100-500/am SEDs of KINGFISH galaxies, Mj(MBB), and those obtained using 
a dust emission model, M d (model). In the left panels, M d (model) is obtained by D12 applying the DL07 method to the whole dust SED. In the 
central panels, M d (model) comes from a simplified application of the DL07 approach to just the Herschel 100-500//m SED; the same in the right 
panels, but using the dust properties from Cll (see text for details). In the top panels, M d (MBB) is derived using a power-law fit to the dust 
absorption cross section from D03 (left and center) and from Cll (right). In the bottom panels, the spectral index is changed from the fitted value 
to/3 = 1.5. Dashed lines show the median value of the ratios. Open circles (and the dotted line) refer to the case for p = 1.5 after correcting the 
/c abs normalization (see Sect. 3 for details). 



quantity sometimes referred to as emissivity). A power law was 
assumed for k&s, 

Kabs = Kabs(^o) X I — I . (2) 

Fits were performed using the procedure of Magrini et al. (20 11): 
since the flux densities of D 12 are not color-corrected, the MBB 
was integrated over the the PACS and SPIRE filter response 
functions before comparing it to the data; for SPIRE, the filter 
response functions for extended sources was selected. The MP- 
FIT IDL^f 2 minimization routines were used (Markwardt 2009). 
Recently, concern has been raised on the use of standard ;f 2 min- 
imizations techniques when fitting MBBs to observed SED, with 
biases arising when both r d and /? are derived, resulting in a spu- 
rious Td-/3 anticorrelation that might (or might not) conceal real 
grain properties (see e.g. Kelly et al. 2012). However, in this 
work I always keep f3 fixed and only derive and r d , in anal- 
ogy to the other works I compare with. Results from the fit where 
checked with a bootstrapping technique, by fitting one hundred 
random representations of each SED (each compatible with the 
original photometry, within its error). The mean and standard 
deviations for M d and r d obtained with the bootstrapping were 
almost the same as those obtained by fitting the observed SED 
with MPFIT. For the photometric errors provided by D12, the 
relative errors obtained by the fit on Md and T d are, on average, 
8 and 2%, respectively. 

Since the DL07 emission model used by D12 is based on 
Draine's (2003; hereafter, D03) model for Milky Way (MW) 
dust , an appropriate choice for k r \, s is the absorption cross sec- 
tion of the latter, averaged over the the grain size distribution and 
composition. The Ry = 3.1 MW dust model 1 absorption cross 

1 The averaged absorption cross section for this model is available 
at: http://www.astro.princeton.edu/ draine/dust/dustmix.html. Further 



section for 70/ira < A < IQQjim is well fitted by Eq. 2 with 
/f abs (250//m) = 4.0 cm 2 g -1 , /? = 2.08. (3) 

The ratio between the dust masses obtained with the two 
methods is shown in Fig. 1 (upper-left panel). As in the anal- 
ogous Fig. 9 of D12, the ratio is plotted versus the 70//m-to- 
160//m flux density ratio. The dust mass obtained by fitting a 
MBB to the 100/im < A < 500/im data is quite close to that 
obtained with the DL07 model. The median ratio is 0.96 2 . 

D12 claimed that the dust masses derived with a j3 = 2 
MBB fit are underestimated, on average, by 25% (though it ap- 
pears that they have included the 70/im flux density in the fits, 
which I don't use here). However, they adopted A- a b s (250//m) = 
4.8 cm 2 g -1 , a value 20% higher than that of the more recent 
D03's model (i.e. they underestimated the MBB masses by 
20%). If they had used a cross section consistent with the dust 
model used within their implementation of the DL07 method, 
their result would have been in line with what I find here. 

updates to the model parameters (DL07; Aniano et al. 2012) did not 
change substantially the mean dust absorption cross section for the 
wavelength range considered here (Draine, private communication). 
2 In a number of Herschel papers (see e.g. Magrini et al. 201 1; Davies 
et al. 2012; Smith et al. 2012), it was used 

K abs (350yum) = 1.92 cm 2 g _1 , yS = 2, 

where the normalization at 350yimi is the value for the = 3.1 Milky 
Way (MW) averaged dust absorption cross section from D03. The me- 
dian ratio obtained using these values is 0.95, almost the same of that 
obtained using the fit of Eq. 3. When the exact, tabulated, averaged 
cross section from D03's model is used, the ratio becomes 0.98. It is to 
be noted that I consider here the full dust mass obtained by D12. If, in 
addition to the use of the tabulated values, I correct for the fraction of 
dust that is heated by higher intensity radiation fields (as given by the y 
parameter provided by D12) and does not contribute significantly to the 
emission at A > 100//m, the median ratio rises to 1.0. 
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The dispersion of the mass ratios is 0.18. As it can be seen 
from Fig. 1, this is higher than the error on the ratio, whose 
mean is 0.12 (estimated using the error from the fits for indi- 
vidual MBB masses, and assuming the same errors for the DL07 
masses, which were not provided). The large scatter is related 
to the intrinsic differences of the two methods: the masses from 
D12 were derived using additional fluxes, not just the five Her- 
schel datapoints, with emission from dust heated by the mean 
ISRF constrained by shorter wavelength observations also; the 
dust absorption cross section also depend, though slightly for 
the Herschel wavelength range, on the fraction of the dust mass 
in the form of polycyclic aromatic hydrocarbons (PAHs), the pa- 
rameter q PAli which is a result of D12 fits. 

Indeed if these differences are removed, the scatter is re- 
duced. This was tested deriving DL07 masses using a simplified 
approach: I have used the emission models 3 for dust heated by a 
single intensity of the ISRF only (as defined by the U parameter 
in DL07, a scale factor with respect to the local Galactic ISRF 
from Mathis et al. 1983); and those computed for the Ry = 3.1 
MW dust model only (i.e. q PAil = 4.68%; DL07). As in the 
full DL07 approach, emission templates for different intensites 
were convolved with the filter response functions, but only for 
the five Herschel datapoints considered for the MBB fits. SED 
fits were produced by minimizing x 2 over the provided U grid 
and over and errors on the two quantities were computed 
with the same bootstrapping technique used for MBBs. This 
confirmed that the errors on are similar for both the MBB and 
DL07 approach, as assumed earlier in this section. The results 
of the fits are compatible with those of the full analysis in D12, 
though the Us found here are generally lower than their equiv- 
alent [/nun (a distribution of Us is used in the full approach, the 
dust heated by U m i n being responsible for most of the FIR peak). 
As a consequence, the DL07 masses derived here are larger, but 
only by ~ 10%, the median ratio is smaller, 0.9, and the scatter 
is reduced to 0.07, of the same order of the error (Fig. 1, upper- 
middle panel; the ratio is computed using the same MBB masses 
of top-left panel). This result is identical to the tests we did in 
Magrini et al. (2011). 

The DL07 approach can be used with any dust grain model. 
As a further test, I use it here with the MW dust model by Com- 
piegne et al. (201 1, hereafter, CI 1), which basically differs from 
D03's model in the use of optical properties for amorphous car- 
bon. Using the DustEM code described in that paper 4 I have 
computed the emission for a grid of U values, and derived the 
dust masses as described in the previous paragraph. DustEM 
also provides the mean absorption cross section of the dust dis- 
tribution, which can be fitted with Eq. 2 and 

K abs (250,«m) = 5.1 cm 2 g" 1 , J3 = 1.91. (4) 

As a results of the different K aDS , dust masses derived using the 
simplified, DL07 approach with the CI 1 model are a factor 0.73 
the corresponding masses using D03. When MBB masses are 
computed using Eq. 4, the median ratio between MBB and DL07 
masses for the CI 1 model is 0.82, with a scatter of 0.05 5 (Fig. 1, 
upper-right panel). 

Thus, provided that an appropriate /e a b s is used, MBB fits can 
retrieve the dust mass within at most 20% of what more complex 
models including a distribution of dust grains (and temperatures) 
can. This difference is smaller than the current uncertainties in 

3 The original DL07 models are available at: 
http://www.astro.princeton.edu/ draine/dust/irem.html. 

4 Available at: http://www.ias.u-psud.fr/DUSTEM/ 

5 It raises to 0.84 if the tabulated absorption cross section is used. 



MW A^bs models themselves. Finally, it is worth noting that for 
some objects SEDs are well fitted by neither a MBB, nor a DL07 
model (D12). 

3. MBB dust masses vs fi 

D12 also made MBB fits allowing /? to vary. For most of the 
objects, they find /? ~ 1.5. The analysis of Herschel colors of 
galaxies does suggest that 1^/3^2 (see e.g Boselli et al. 2012; 
Auld et al. 2012). Keeping the same cross section normalization, 
D12 find that dust masses can be severely underestimated. In- 
deed, if I use the normalization of Eq. 3 together with /? = 1 .5 in 
Eq. 2, the median ratio between MBB and DL07 masses reduces 
to 0.55, 0.51 and 0.54 for the three comparison in Fig. 1 (see the 
lower panels, filled symbols). This is due to the change in the fit- 
ted r d , whose median value raises from »20K for /? = 2.08 and 
21K for/? = 1.91 to «24 for/? =1.5 (for a fixed normalization of 
/c a bs, a smaller amount of hotter dust can reproduce the observed 
fluxes). 

Despite it being common practice, it is however not justi- 
fyied to use /c a bs(^o) from a dust model (which has a proper 
spectral index, close to /? = 2 for both models considered here) 
and then assume a different spectral index. This inconsistency 
leads to the puzzling result that the dust mass estimate can vary 
with the wavelength chosen for the normalization: for the case 
of /? — 1.5, adopting the model cross section at 500//m results 
in dust masses that are 50% to 30% higher than those obtained 
with the 250//m normalization of Eq. 3 and 4, respectively. The 
difference simply comes from (500/250)^ L5) , the fitted temper- 
ature being independent of the choice for the normalization. A 
similar remark was made by Skibba et al. (2011), which fitted 
KINGFISH SEDs from Spitzer and SPIRE using /? = 1.5 and 
/e a b s (500yum) from D03. However, they inexplicably concluded 
that their normalization results in a mass that is a factor 3 higher 
than when /e a b s (250/im) is used. The temperatures from Skibba 
et al. (2011) are sistematically higher than those obtained here 
for/? = 1.5, and the dust masses, when corrected to the distances 
and normalization used here, are lower by a factor 2.7; A similar 
discrepancy was also noted by Galametz et al. (2012) and im- 
puted to the use of 70//m Spitzer data, which is contaminated by 
emission from stochastically heated grains. 

For a proper mass estimate, one would need an absorption 
cross section derived consistently with the assumption made on 
/?. This could come from a galaxy whose SED behaves as a 
MBB with the chosen /?, provided the dust mass can be derived 
with sufficient accuracy independently of the SED itself (for ex- 
ample using the mass of metals as a proxy, as done by James 
et al. 2002); or from a dust model for the MW grains in which 
the materials cross sections follow the chosen spectral behav- 
ior (provided that those materials are able to predict consistently 
both the MW extinction curve and the FIR/submm emission). 

As a numerical experiment, I show here what a different 
choice of /? could imply in the derivation of the absorption cross 
section from the MW SED. Cll derived the FIR/submm sur- 
face brightness for the Diffuse High Galactic Latitude medium 
(DHGL, with latitude \b\ > 15°) normalised to the hydrogen 
column density, I v /Nn- They provide Herschel surface bright- 
nesses for the DHGL by convolving their COBE-FIRAS spec- 
trum with the PACS (160/im) and SPIRE (extended emission) 
filter response functions. They also provide the COBE-DIRBE 
100/luti surface brigthness, which I converted to PACS (lOOyum) 
by multiplying it by 1.1 to take into account the different filter 
widths. This dataset is thus analogous to what I used in Sect. 2 
to derive dust masses. 
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Fig. 2. MBB fits to the DHGL emission. 

For a dust grain of temperature emitting as a MBB with 
a power-law absorption cross section (i.e. the same assumptions 
as in Sect. 2), the normalization of the cross section can be de- 
rived by fitting 



h Tabs D v , 
where 



d >' 



Tabs rr>lf-\ T. dbs (Ao) 
— = Kabs m H (D/G)mw = — X | 



(6) 



tMW 
1 d 



The absorption cross section per H atom, t^s/Nh, can be con- 
verted into the cross section per unit dust mass of Eq. 2 by divid- 
ing it for the hydrogen nucleon mass m H and the MW dust-to-gas 
mass ratio (D/G)mw- 

Using /3 = 1.91, the DHGL emission can be fitted 6 with 
= 17.8 ± 0.4K. At 250/im, the absorption cross section 
is found to be T. dbs (250fim)/N H = 0.84 + 0.09 x 10~ 25 cm 2 H _1 . 
Assuming (D/G)mw ~ 0.01, as can be estimated from the ele- 
mental depletion patterns in the diffuse MW gas (Draine 201 1), 
it is K a bs(250//;7i) ss 5.0 cm 2 g _1 , a value close to that from 
Cll in Eq. 4: this is not surprising, since the dust model was 
made to fit the DHGL data and has a similar dust-to-gas ra- 
tio, (D/G)mw = 0.0102. For = 2, T™ w = 17.4 + 0.4K 
and Tab s (250 y um)/A^ H = 0.91 ± 0.10 x 10~ 25 cm 2 H 1 (Fig. 2), 
a value close to the original determination of Boulanger et al. 
(1996) on the COBE-F1RAS spectrum (though Cll correct for 
the contribution of ionised gas to the hydrogen column den- 
sity, while Boulanger et al. only consider atomic gas). When 
(3 = 2.08, as for D03's model, Tf* = 17.1 ± 0.4K and 

Tabs(250//m)/jV H = 0.98 + 0.11 x 10~ 25 cm 2 H _1 , which con- 
verts to K a bs(250yum) x 5.9 cm 2 g _1 , a value 45% higher than 
that in Eq. 3. The model could still be reconciled with the data, 
if the ionised gas contributes more to the total hydrogen column 
density than what assumed by Cll. To be consistent with the 
mass determinations done so far, I continue to use the absorp- 
tion cross sections from Eqs. 3 and 4 but scale them according to 
the ratios of Tabs/A^ obtained from the DHGL with different /3s. 



6 I stress here that the aim of this numerical tests is only to show how 
the results are affected by the choice of /?, and not to derive the spec- 
tral index of MW emission. Though /?s between 1.5 and 2.1 produce 
plausible fits in the limited spectral range considered in this work, the 
analysis of the full MW dust SED yields « 1.8 ± 0.2 (Planck Collab- 
oration et al. 2011). 



Using = 1.5 the temperature rises to Tf w = 19.7 + 0.5K 
and the cross section reduces to T a bs(250yum)/7V H = 0.57 + 0.06 x 
10~ 25 cm 2 H _1 , a factor 0.58 and 0.68 those for (3 = 2.08 and 
1.91, respectively. Obviously, the change of ft has the same ef- 
fect on both the dust mass and cross section determination, since 
Eq. 1 and 5 are formally identical. 

When the normalizations of Eq. 3 and 4 are corrected by 
those factors, the dust masses obtained using (3 = 1 .5 become 
very close to those obtained using (3 fitted to the dust models 
(see the open symbols in the lower panels of Fig. 1, with median 
mass ratios of 0.94, 0.88 and 0.8 for the three panels, from left to 
right). Thus, as already noted in Bianchi et al. (1999), the dust 
mass estimate does not change significantly with (3 if both the 
absorption cross section and the dust mass are derived with the 
same f3. Substituting Eq. 2, 5 and 6 in Eq. 1, one finds that the 
dust mass of a galaxy depends on f3 only through the ratio 



M d 



B v (r d MW ) 

B V (T A ) ■ 



For the wavelength and temperature ranges considered here, this 
ratio does not change significantly if the temperatures are de- 
rived using j3 = 1 .5 or j3 x 2. 



(5) 4. Conclusions 



I have shown that, if the DL07 model of dust emission in galax- 
ies is correct, a simple single temperature MBB fit to the SED for 
A > lOOyum can provide in a reliable way one of the parameters 
of the more complex approach, i.e. the dust mass. For this, it is 
necessary to ensure the consistency in dust emission properties 
between the two approaches: the absorption cross section used in 
the MBB needs to be the average over the grain size distribution 
and composition of the dust model used in DL07. This consis- 
tency is broken when a cross section normalization at a given 
wavelength is taken from, e.g., the D03 and Cll grain models 
(which have been derived with /3 ~ 2, implicity coming from the 
adopted material properties) and used in conjunction with /3 & 2 
in the MBB fits. 

The need for consistency is illustrated by a numerical exper- 
iment: using either (3 » 2 or 1.5 to derive the absorption cross 
section from the MW SED for 100/j.m < A < 500yum, and the 
same ft to derive the dust mass from the SED of galaxies (both 
under the assumption of a single temperature MBB) results in 
dust masses which are almost independent of /3. This might be 
true only for this limited range, which is still able to provide a 
reliable fit of the MW SED. For /3s outside this range, the nor- 
malization of the absorption cross section must come from other 
sources. 

Finally, I remind the reader that I have only considered the 
SED for A > 100/Lmi. Single temperature MBB fits to datasets 
including flux densities at shorter wavelengths (i.e. with a large 
contribution of non-equilibrium emission), as well as two tem- 
peratures MBB fits, might result in further biases in the deriva- 
tion of the dust mass, in addition to those discussed here (Aniano 
et al. 2012, Aniano et al. in preparation). 

Acknowledgements. I thank Daniel Dale for providing me the parameters of the 
DL07 fits, Bruce Draine for stimulating discussions, Carlo Giovanardi and Leslie 
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